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—  A  one-dimensional  model  (CE-QUAL-Rl)  of  water  quality  for  reservoirs  was 
evaluated  using  data  collected  on  the  Eau  Calle  Reservoir,  in  west-central 
Wisconsin.  Data  collected  in  1981  were  used  for  calibration;  data  collected 
in  1982  were  used  for  confirmation.  Graphical  and  statistical  comparisons 
were  made  for  over  6,500  samples  representing  21  variables.  The  statistical 
test  used  for  comparing  measured  versus  predicted  values  was  the  Reliability 
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PREFACE 


This  report  was  sponsored  by  the  Office,  Chief  of  Engineers  (OCE) , 
US  Army,  as  part  of  the  Environmental  and  Water  Quality  Operational 
Studies  (EWQOS),  Work  Units  IB.l  and  IC.l.  OCE  Technical  Monitors  for 
EWQOS  were  Dr.  John  Bushman,  Mr.  Earl  Eiker,  and  Mr.  James  L.  Gottesman. 

Confirmation  of  CE-QUAL-R1  using  data  from  Eau  Galle  Reservoir  and 
the  writing  of  this  report  were  accomplished  by  Dr.  Joseph  H.  Wlosinski 
and  Dr.  Carol  D.  Collins,  Water  Quality  Modeling  Group  (WQMG) ,  Ecosystem 
Research  and  Simulation  Division  (ERSD) ,  of  the  Environmental  Laboratory 
(EL),  US  Army  Engineer  Waterways  Experiment  Station  (WES).  Total  con¬ 
firmation,  other  model  development,  and  data  collection  and  manipulation 
were  performed  by  many  other  members  of  the  ERSD  during  the  course  of 
EWQOS.  The  draft  report  was  reviewed  by  Dr.  James  L.  Martin  and 
Dr.  Stephen  P.  Schreiner,  both  of  the  WQMG.  The  report  was  edited  by 
Ms.  Jessica  S.  Ruff  of  the  WES  Publications  and  Graphic  Arts  Division. 

The  study  was  conducted  under  the  direct  supervision  of 
Mr.  Mark  S.  Dortch,  Chief,  WQMG,  and  under  the  general  supervision  of 
Mr.  Donald  L.  Robey,  Chief,  ERSD,  and  Dr.  John  Harrison,  Chief,  EL,  WES. 
Program  Manager  of  EWQOS  was  Dr.  Jerome  L.  Mahloch,  EL. 

During  the  study  and  preparation  of  this  report,  COL  Tilford  C. 
Creel,  CE,  and  COL  Robert  C.  Lee,  CE,  were  Commanders  and  Directors  of 
WES  and  Mr.  Fred  R.  Brown  was  Technical  Director.  At  the  time  of  publi¬ 
cation,  COL  Allen  F.  Grum,  USA,  was  Director  and  Dr.  Robert  W.  Whalin 
was  Technical  Director. 

This  report  should  be  cited  as  follows: 

Wlosinski,  J.  H. ,  and  Collins,  C.  D.  1985.  "Confirma¬ 
tion  of  the  Water  Quality  Model  CE-QUAL-R1  Using  Data 
from  Eau  Galle  Reservoir,  Wisconsin,"  Technical  Report 
E-85-11,  US  Army  Engineer  Waterways  Experiment  Station, 

Vicksburg,  Miss. 


CONFIRMATION  OF  THE  WATER  QUALITY  MODEL  CE-QUAL-R1 
USING  DATA  FROM  EAU  GALLF.  RESERVOIR,  WISCONSIN 


PART  I:  INTRODUCTION 


Background 


1  .  CE-QUAL-R1  is  a  water  quality  model  that  is  undergoing  contin¬ 
uing  development  by  the  Corps  of  Engineers  as  part  of  the  Environmental 
and  Water  Quality  Operational  Studies  (EWQOS).  The  model  predicts  the 
vertical  distribution  of  temperature  and  other  water  quality  variables 
in  a  reservoir  through  time.  The  version  of  the  model  used  for  this 
study  was  basically  the  one  described  in  the  revised  User's  Manual  (En¬ 
vironmental  Laboratory  1982),  with  changes  as  reported  in  Wlosinski  and 
Collins  1985.  In  addition,  a  weir  subroutine  and  a  macrophyte  subrou¬ 
tine  were  added  to  the  model  to  better  simulate  conditions  of  the  Eau 
Galle  Reservoir.  All  work  was  done  on  a  VAX-11-750  computer. 

2.  Model  confirmation  using  data  collected  on  the  Eau  Galle  Res¬ 
ervoir  follows  a  similar  study  using  data  collected  on  DeGray  Lake,  a 
Corps  of  Engineers  project  located  in  the  Ouachita  Mountains  in  south- 
central  Arkansas  (Wlosinski  and  Collins  1985).  One  purpose  of  the  two 
studies  was  to  test  the  model  on  reservoirs  that  differed  markedly. 

DeGray  Lake  is  32  km  long,  dendritic,  and  has  a  normal  depth  near  57  m, 

7  2  8  3 

an  area  of  5.3  *  10  m  ,  and  a  volume  of  7.9  *  10  m  .  Eau  Galle  Reser¬ 
voir  is  1  km  long,  circular,  and  has  a  normal  depth  near  10  m,  an  area 
5  2  6  3 

of  4.5  x  10  m  ,  and  a  volume  of  1.7  k  10  m  .  DeGray  is  monomictic, 
while  Eau  Galle  is  dimictic.  Biologically,  Eau  Galle  is  much  more  pro¬ 
ductive  than  DeGray. 


Purpose 

3.  The  purpose  of  this  study  was  to  ensure  that  model  performance 
is  suitable  for  the  needs  of  Corps  of  Engineers  (CE)  District  and 


Division  Offices.  This  was  accomplished  by  comparing  model  predictions 
to  field  measured  values.  Data  collected  in  1981  were  used  for  model 
calibration,  while  data  collected  in  1982  were  used  for  confirmation. 
Graphical  and  statistical  tests  were  used  for  the  comparisons,  which 
consisted  of  the  concentrations  of  state  variables.  In  addition,  the 
flux  between  modeled  components  was  compared  in  cases  where  data  were 
available . 


Model  Description 

4.  CE-QUAL-R1  is  a  water  quality  model  that  computes  the  vertical 
distribution  of  37  physical,  biological,  and  chemical  constituents  in  a 
reservoir  through  time.  In  the  model,  a  reservoir  is  conceptualized  as 
a  vertical  series  of  horizontal  layers  where  thermal  energy  and  mass  are 
uniformly  distributed  in  each  layer.  The  horizontal  layer  thicknesses 
are  variable  and  dependent  on  the  balance  of  inflowing  and  outflowing 
waters.  Variable  layer  thicknesses  permit  accurate  mass  balancing  and 
reduce  numerical  dispersion  during  periods  of  large  inflow  and  outflow. 

5.  Inflowing  waters  are  distributed  vertically  based  on  density 
differences,  so  that  simulations  of  surface  flows,  interflows,  and 
underflows  are  possible.  Water  density  is  dependent  on  temperature  and 
dissolved  and  suspended  solids  concentrations.  Outflowing  waters  are 
withdrawn  from  the  horizontal  layers  considering  density  stratification 
using  the  selective  withdrawal  algorithm  of  Bohan  and  Grace  (1973). 
Reservoir  outflows,  by  port,  can  either  be  specified  or  the  user  can 
invoke  a  subroutine  that  will  choose  port  flows  in  order  to  meet  a  down¬ 
stream  temperature  objective. 

6.  The  heat  budget  includes  the  components  of  short-  and  long¬ 
wave  radiation,  back  radiation,  reflected  solar  and  atmospheric  radia¬ 
tion,  evaporative  loss,  conductive  heat  transfers,  and  gain  or  loss 
through  inflows  and  outflows.  Vertical  transport  of  thermal  energy  and 
mass  is  achieved  through  entrainment  and  turbulent  diffusion.  Entrain¬ 
ment  determines  the  depth  of  the  upper  mixed  layer  and  the  onset  of 
stratification.  It  is  calculated  from  the  turbulent  kinetic  energy 


influx  generated  by  wind  shear  and  convective  mixing  using  an  integral 
energy  approach  (Johnson  and  Ford  1981).  Turbulent  diffusion  tends  to 
reduce  gradients  and  is  incorporated  using  a  turbulent  or  eddy  diffusion 
coefficient  that  is  dependent  on  the  windspeed,  magnitude  of  inflows  and 
outflows,  and  density  stratification. 

7.  The  prediction  of  water  quality  is  based  upon  simulation  of 
the  interaction  of  numerous  biological  and  chemical  constituents. 

Forces  that  directly  affect  the  simulation  of  these  constituents  are 
temperature,  irradiation,  wind  speed,  inflow  and  outflow  rates,  and 
inflowing  and  outflowing  masses.  The  physical  distribution  of  mass  is 
dependent  upon  the  diffusive  and  convective  processes  described  above 
and  on  settling  processes. 

8.  The  interactions  between  and  among  the  constituents  ulti¬ 
mately  determine  the  constituent  masses.  Photosynthesis,  dark  respira¬ 
tion,  photorespiration,  and  nonpredatory  mortality  are  processes  influ¬ 
encing  algal  and  macrophyte  mass.  Grazing  by  fish  and  zooplankton  are 
additional  influences  on  algae.  Ingestion,  egestion,  and  respiration 
affect  zooplankton  and  fish  growth.  Inorganic  compounds  such  as 
ammonia-N,  nitrite  plus  nitrate-N,  orthophosphate,  and  silica  are  con¬ 
sumed  and  produced  as  a  result  of  the  photosynthetic  and  respiratory 
processes  of  the  plants  and  animals.  Orthophosphate  and  ammonia-N  are 
adsorbed  to  solids  according  to  a  modified  equation  for  the  Langmuir 
isotherm.  Ammonia-N  is  also  removed  by  conversion  to  nitrite  plus 
nitrate-N  under  aerobic  conditions.  Nitrite  plus  nitrate-N  is  lost 
through  denitrification. 

9.  Hass  of  detritus  is  dependent  on  algal  and  macrophyte  mortal¬ 
ity,  ingestion  by  fish  and  zooplankton,  egestion  of  zooplankton,  and 
settling.  Decomposition  of  detritus  contributes  mass  to  ammonia-N, 
nitrite  plus  nitrate-N,  orthophosphate,  and  inorganic  carbon. 

10.  Inflowing  and  initial  concentrations  for  dissolved  organic 
matter  (DOM)  are  fractionized  between  labile  and  refractory  DOM  compart¬ 
ments.  Refractory  DOM  is  usually  more  resistant  to  decomposition  than 
labile  DOM  and  decomposes  at  a  slower  rate.  Plant  photorespiration  con¬ 
tributes  to  labile  DOM.  As  labile  DOM  decomposes,  products  are 
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distributed  to  inorganic  nutrients  and  refractory  DOM. 

11.  Dissolved  oxygen  concentration  is  of  primary  importance  to 
reservoir  management.  Oxygen  is  evolved  by  algal  and  macrophyte  photo¬ 
synthesis.  Oxygen  demand  in  CE-QUAL-R1  is  created  by  the  processes  of 
nitrification,  decomposition  of  organic  compounds  and  sediment,  respira¬ 
tion,  and  oxidation  of  reduced  products  of  anaerobic  reactions.  Oxygen 
may  also  be  gained  or  lost  at  the  air-water  interface.  Anaerobic  and 
aerobic  conditions  resulting  from  changes  in  oxygen  concentration  drive 
many  other  modeled  processes.  If  the  system  becomes  anaerobic,  decompo¬ 
sition  of  organic  material  slows  considerably,  and  certain  compounds  are 

released  from  the  sediment.  These  compounds  include  ammonia-N 
+  -3 

(NH4  -N) ,  orthophosphate-P  (P04  -P),  dissolved  reduced  manganese  and 

+2  +2  -2 

iron  (Mn  ,  Fe  ),  and  sulfide  (S  ).  S'  ments  release  almost  all  the 
anaerobic  compounds  generated  in  CE-QUAL-R1 ;  reductions  and  inflow 
account  for  the  remainder.  Reoxygenation  of  the  system  will  reverse 
these  reactions. 

12.  Total  dissolved  solids  are  simulated  to  obtain  an  approxima¬ 
tion  of  ionic  strength.  Calculations  based  on  the  equilibrium  reactions 
of  bicarbonate,  carbonate,  and  hydroxyl  ions  and  on  ionic  strength 
result  in  the  pH  value  reported  for  each  layer.  This  value  is  then  used 
to  calculate  the  carbon  dioxide  concentration  which  contributes  to  plant 
growth  and  diffuses  across  the  air-water  interface.  Total  alkalinity  is 
simulated  in  CE-QUAL-R1  to  provide  an  indication  of  the  buffering  capac¬ 
ity  of  the  system.  Alkalinity  is  modeled  as  a  conservative  substance, 
being  only  advected  and  diffused. 

13.  Suspended  solids  influence  both  the  density  and  light 
regimes.  Suspended  solids  are  subjected  to  advection,  diffusion,  and 
settling.  A  more  detailed  description  of  the  final  model  used  in  this 
study  will  be  available  in  a  revised  user’s  manual. 

Eau  Galle  Reservoir  Study  Site 

14.  Eau  Galle  Reservoir  is  located  in  west-central  Wisconsin, 

56  km  upstream  of  the  confluence  of  the  Eau  Galle  and  Chippewa  Rivers 
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and  immediately  upstream  of  the  town  of  Spring  Valley.  The  drainage 

2 

area  at  the  damsite  is  165  km“,  with  the  majority  of  the  land  being  used 
for  dairy  operations  and  associated  agriculture.  At  a  normal  pool  ele¬ 
vation  of  286.5  m  mean  sea  level  (1929  adjustment),  the  surface  area  of 
2  6  3 

the  pool  is  0.5  km  ,  the  volume  is  1.7  x  10  m  ,  the  maximum  depth  is 
9.75  m,  and  the  reservoir  mean  depth  is  3.7  m.  The  dam  is  a  rolled 
earth-fill  structure,  with  a  length  of  486  m  and  a  maximum  height  above 
the  streambed  of  37  m.  Pool  elevation  is  controlled  by  a  vertical  slide 
gate  leading  to  a  horseshoe  conduit  with  an  upstream  invert  elevation  of 
279.3  m  and  by  an  uncontrolled  morning  glory  weir  at  286.5  m.  During 
1982  and  1983,  the  pool  fluctuated  between  286.3  and  288.6  m  although  it 
rarely  exceeded  286.8  m.  A  more  detailed  description  of  the  dam,  reser¬ 
voir,  and  watershed  is  available  in  Ashby  (1985). 

15.  An  overview  of  the  reservoir  is  presented  in  Figure  1.  The 
dam  is  located  in  the  top  right  corner  of  the  figure.  The  in-pool 
sampling  stations  are  labeled  10,  20,  30,  40,  50,  and  60.  The  Eau  Galle 
River  and  French  Creek  join  just  upstream  of  the  reservoir,  and  enter  at 
the  point  labeled  A.  Lousy  and  Lohn  Creeks  enter  at  B  and  C, 
respectively. 
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PART  II:  DATA  REQUIREMENTS 

16.  Four  types  of  data  were  required  for  the  simulation  studies: 
initial  conditions,  driving  variables  (also  termed  boundary  conditions 
or  updates),  model  coefficients,  and  conf irmation  data  (also  termed  cal¬ 
ibration  or  verification  data).  Initial  conditions  refer  to  the  amounts 
of  state  variables  at  the  start  of  simulation.  Driving  variables  act  to 
change  the  state  variables  of  the  system,  but  lie  outside  the  system. 

For  CE-Ql'AI.-R  1 ,  the  driving  variables  include  meteorological  data,  flow 
from  upstream,  and  controlled  outflow  releases.  Model  coefficients  are 
constants  used  in  the  algorithms  that  comprise  the  model.  Confirmation 
data  refer  to  field  measured  values  which  are  then  compared  to  model 
predictions.  Complete  data  sets  are  included  in  Appendix  A  for  1981  and 
for  1981.  These  data  sets  include  only  those  data  needed  to  simulate 
the  Eau  dalle  Reservoir.  Confirmation  data  are  not  included.  Informa¬ 
tion  on  the  sampling  program  and  analytical  methods  can  be  found  in 
Johnson  and  l.auer  (1985). 


Initial  Conditions 


17.  Initial  conditions  were  measured  on  7  April  1981  and  20  April 
1982  at  the  deepest  part  of  the  pool,  located  at  station  20  (Figure  1). 
Those  variables  for  which  values  were  available  were  temperature,  alka¬ 
linity,  ammonia-N,  nitrite  plus  nitrate-N,  oxygen,  orthophosphate-P, 
total  dissolved  solids,  pH,  suspended  solids,  silica,  coliforms,  sul¬ 
fate,  and  sulfide.  Reduced  and  oxidized  manganese  and  iron  were  ob¬ 
tained  from  dissolved,  and  from  the  difference  between  total  and  dis¬ 
solved,  manganese  and  iron,  respectively. 

18.  Initial  conditions  for  algae  were  estimated  using 
chlorophy 1 1 -a  measurements  end  information  from  Barko  et  al .  (1985). 

They  found  a  significant  relationship  (r2  =  0.84)  of  2.8  mg 

ch 1 orophv 1 1-a  per  gram  of  fresh  weight.  Using  this  value,  and  the 
assumption  that  90  percent  of  fresh  weight  is  water,  the  conversion  to 


modeled  units  of  grams  dry  weight  was  made.  Barko  et  al.  (1985)  also 
found  that  algal  populations  were  dominated  by  Bacillariophyceae  (dia¬ 
toms),  Cyanophyta  (blue-greens),  and  Pyrrophyta  (dinof lagellates)  ; 
therefore,  this  grouping  was  used  for  the  three  modeled  algal  compart¬ 
ments.  Because  diatoms  dominated  the  reservoir  in  the  spring,  the  total 
chlorophyll  value  for  initial  conditions  was  assumed  to  be  diatoms.  For 
purposes  of  confirmation,  the  mass  of  the  three  algal  compartments  was 
summed  and  converted  to  chlorophvll-a. 

19.  Dissolved  organic  carbon  was  used  to  estimate  refractory  and 
labile  DOM,  with  30  percent  estimated  to  be  labile.  Dissolved  organic 
carbon  was  assumed  to  be  46  percent  of  DOM.  Estimates  of  fish  mass  were 
obtained  from  Leidy  and  Jenkins  (1977).  Because  the  water  was  well 
aerated  on  the  day  that  initial  conditions  were  measured,  values  for 
iron  sulfide  were  assumed  to  be  zero.  The  other  variables  that  were  not 
measured,  detritus,  zooplankton,  and  materials  in  the  sediments,  were 
given  low  initial  values  or  values  based  on  the  DeGray  Lake  study.  The 
ranges  of  values  for  initial  conditions  for  the  two  simulation  years  are 
listed  in  Appendix  B. 


Driving  Variables 

20.  Initial  simulations  used  meteorological  data  measured  at  the 
National  Weather  Service  station  at  Minneapolis-St .  Paul,  Minn.,  80  km 
west  of  the  Eau  Galle  Reservoir.  Daily  averages  were  used  for  dry  bulb 
and  dew  point  temperatures,  cloud  cover,  barometric  pressure,  and  wind 
speed.  During  calibration  of  the  1981  data  set,  it  was  noticed  that 
there  was  too  much  mixing  during  a  few  critical  periods  of  the  year.  A 
comparison  of  wind  speeds  between  the  Minneapolis-St.  Paul,  Minn.,  and 
Eau  Claire,  Wis.,  weather  stations  showed  a  fairly  consistent  difference 
between  the  stations,  with  the  Eau  Claire  values  being  lower.  The  Eau 
Claire  station  is  approximately  68  km  east  of  the  reservoir.  For  exam¬ 
ple,  for  the  first  week  of  May  1981,  wind  speeds  averaged  17.8  km/hr  at 
the  Minneapolis  station,  and  13.5  km/hr  at  Eau  Claire,  with  maximum 
winds  for  the  period  of  31.3  and  23.2  km/hr,  respectively.  The 


calibration  simulation  was  much  better  with  the  lower  wind  values  from 
Eau  Claire.  The  1982  confirmation  simulation  showed  the  same  pattern; 
therefore,  wind  data  from  St.  Paul  were  replaced  with  values  measured  at 
Eau  Claire.  The  wind  data  used  for  the  simulations  are  shown  in  Fig¬ 
ure  2  (1981)  and  Figure  3  (1982).  The  rest  of  the  meteorological  data 
came  from  the  St.  Paul  station. 
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Figure  2.  Wind  speed  used  for  the  1981  Eau  Galle  data  set 

21.  Two  tributaries  were  used  for  simulation  purposes.  One  tribu¬ 
tary  represented  the  Eau  Galle  River  and  French  Creek;  the  other  repre¬ 
sented  Lousy  and  Lohn  Creeks.  For  the  first  simulation,  inflow  and  out¬ 
flow  discharge  values  were  taken  from  US  Geological  Survey  (USGS) 
records  (1981,  1982),  corrected  for  watershed  areas  above  the  gages  when 
compared  with  areas  above  the  reservoir.  These  values  were  not  satis¬ 
factory,  predicting  pool  levels  much  lower  than  measured.  In  checking 
the  USGS  data  for  the  period  November  1980  to  September  1981,  the  addi¬ 
tional  discharge  from  the  reservoir,  corrected  for  areas  above  the 
gages,  was  in  error  enough  to  drain  the  reservoir  four  times  over  at  the 
normal  pool  elevation.  Instead  of  using  the  USGS  data,  inflow  and 
outflow  discharge  values  were  obtained  from  daratender  records,  stored  at 
the  US  Army  Engineer  District  (USAED) ,  St.  Paul.  These  outflow  values 
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JULIAN  DAY 

Figure  3.  Wind  speed  used  for  the  1982  Eau  Galle  data  set 

were  estimated  using  a  tailwater  rating  curve,  and  the  inflow  values 
were  back-calculated  using  the  outflow  values  and  the  change  in  pool 
storage.  To  distribute  the  inflow  between  the  two  modeled  tributaries, 
the  ratio  from  the  USGS  data  was  used. 

22.  Concentrations  of  constituents  in  the  inflowing  water  were 
measured  biweekly  on  all  four  tributaries.  Concentrations  for  the 
modeled  tributaries  were  calculated  based  on  the  measured  value  and  the 
amounts  of  inflow  from  each  tributary.  Orthophosphate  phosphorus  value 
were  found  to  be  related  to  discharge,  so  daily  concentrations  were 
estimated  based  on  flow  using  regression  techniques.  Daily  temperature 
data  were  available  from  USGS  records  measured  at  Woodville,  Wis., 
approximately  5  km  upstream  of  the  reservoir.  Using  these  data,  inflow 
temperatures  measured  for  the  four  tributaries  on  a  biweekly  basis,  and 
regression  techniques,  daily  temperature  estimates  for  the  two  modeled 
tributaries  were  obtained. 


Model  Coefficients 


23.  Coefficients  included  in  physical  and  thermal  equations  of 
the  model  were  calibrated  using  CE-THERM-R1,  the  thermal  portion  of 
CE-QUAL-R1,  and  temperature  data  collected  in  1981.  These  coefficients, 
and  others  originally  estimated  from  literature  values  (Collins  and 
Wlosinski  1983,  Jorgensen  1979)  and  other  model  studies  (Wlosinski  1981, 
Wlosinski  and  Collins  1985)  for  biological  and  chemical  processes,  were 
used  in  CE-QUAL-R1 .  Again,  calibration  of  the  coefficients  was  per¬ 
formed  using  reservoir  data  collected  in  1981.  The  same  coefficients 
were  used  for  the  confirmation  simulation  as  were  used  for  the  final 
calibration  simulation.  Values,  units,  and  an  explanation  of  coeffi¬ 
cients  are  included  as  Appendix  C. 
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Confirmation  Data 


24.  The  majority  of  confirmation  data  was  of  the  same  form  as 
data  for  initial  conditions.  The  data  usually  were  collected  biweekly 
at  1-m  increments  of  depth.  Variables  used  for  confirmation  included 
temperature,  chlorophyll-a,  silica,  total  manganese,  total  organic  car¬ 
bon,  dissolved  organic  carbon,  orthophosphate-phosphorus,  inorganic 
carbon,  ammonia-nitrogen,  nitrite  plus  nitrate  nitrogen,  oxygen,  pH, 
alkalinity,  total  dissolved  solids,  suspended  solids,  total  iron,  sul¬ 
fate,  dissolved  manganese  and  iron,  and  sulfide.  Chlorophyll-a  was  com¬ 
pared  to  the  summation  of  all  three  modeled  algal  compartments.  Pre- 

+2 

dieted  values  for  Mn  and  Mn  (IV)  were  combined  and  compared  to  total 

manganese.  Total  organic  carbon  included  all  three  algal  compartments 

plus  labile  and  refractory  DOM,  detritus,  and  zooplankton.  Dissolved 

organic  carbon  was  compared  to  the  labile  plus  refractory  compartments. 

+2 

Total  iron  included  Fe  (III),  Fe  ,  and  63  percent  of  iron  sulfide. 

Even  though  the  confirmation  simulation  used  data  collected  in  1982, 
similar  comparisons  were  made  in  1981  for  model  calibration. 
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25.  In  addition  to  the  concentrations  of  the  above  variables, 


flux  rates  were  available  for 
production,*  algal  settling,** 
a  laboratory  study  (Gunnison, 


certain  processes.  These  included  algal 
and  sediment  oxygen  demand  obtained  from 
Chen,  and  Brannon  1983). 


*  Personal  Communication,  November  1983,  J.  W.  Barko,  Research 

Biologist,  US  Army  Engineer  Waterways  Experiment  Station,  Vicksburg, 
Miss . 

**  Personal  Communication,  November  1983,  W.  F.  James,  Physical 

Scientist,  US  Army  Engineer  Waterways  Experiment  Station,  Vicksburg, 
Miss . 


■  v  '  t  VH  in V  ~  .1 


«  -V.-  -Jt  _■  .F 


<r  +  -v 


-"  --  r  -r  -v-  <r  . 


PART  III:  METHODS 


Rationale 


26.  Methods  for  model  confirmation  follow  from  the  objective  of 
developing  the  CE-QUAL-R1  model:  to  provide  CE  District  and  Division 
Offices  with  a  means  of  studying  preimpoundment  and  post  impoundment 
water  quality  problems  and  the  effects  of  reservoir  operation  on  water 
quality.  The  goal  of  the  confirmation  process  was  to  supply  the  best 
possible  tool,  within  the  assumptions  specified  for  CE-QUAL-R1,  for  res¬ 
ervoir  water  quality  management.  Two  main  processes  were  used  during 
the  EWQOS  Program  to  evaluate  the  model  (Wlosinski  19SA):  the  first 
tests  the  software  code;  the  second  examines  model  predictions.  Most  of 
the  testing  of  the  software  code  was  performed  in  earlier  studies,  and 
the  results  have  been  reported  by  Wlosinski  and  Collins  (1985).  The 
main  purpose  of  the  Eau  Galle  simulations,  reported  here,  was  to  examine 
model  predictions. 


Variability 

27.  One  of  the  main  assumptions  of  CE-QUAL-R1  is  that  the  reser¬ 
voir  is  spatially  one  dimensional.  The  model  represents  a  reservoir  by 
a  vertical  series  of  layers,  each  of  which  is  horizontally  uniform. 
Several  factors  contribute  to  variability  between  model  and  prototype. 
Near  tributaries,  inflowing  water  affects  measured  concentrations,  and 
the  one-dimensional  assumption  is  weakened.  Macrophyte  beds  in  shallow 
areas  affect  that  part  of  the  system  more  than  deeper,  macrophyte-free 
areas.  In  addition,  the  equations  of  the  model  cannot  be  solved  in 
closed  form;  time  steps  are  needed  to  find  approximate  solutions  to  the 
equations.  Within  each  modeled  time  step  (often  1  day  for  CE-QUAL-R1), 
details  may  be  lost.  The  model  will  predict  one  value  for  each  variable 
in  each  layer  during  a  time  step.  In  the  real  system,  the  actual  range 
of  values  can  be  great.  Consider,  for  example,  Figure  A,  which  is  a 
graph  of  oxygen  values  measured  at  three  different  stations  during  a 
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DISSOLVED  OXYGEN,  gm~3 

Figure  4.  Example  of  oxygen  variability  (shaded 
area)  at  different  locations  during  a  24-hr  period 
at  Eau  Galle  Reservoir.  Data  were  collected  on 
2  September  1981 

2 4  hr  period  on  2  September  1981.  The  values  range  between  17  and 
182  percent  of  saturation  at  the  temperatures  measured.  With  such  high 
variability,  model  evaluation  would  be  very  difficult,  considering  the 
one-dimensional  assumption.  Rather  than  accounting  for  data  variability 
between  stations,  the  single  predicted  value  for  each  variable  at  each 
layer  was  compared  to  the  value  measured  only  at  station  20.  Station  20 
was  chosen  because  it  was  the  deepest  station  and  had  the  most  data. 

The  evaluation  process  would  then  show  whether  or  not  the  model  is  pre¬ 
dicting  the  major  changes  occurring  in  the  reservoir. 

Evaluation  Process 

28.  Both  graphical  and  statistical  tests  were  made  for  predicted 
versus  measured  data.  Qualitatively,  graphs  were  prepared  showing 


vertical  profiles  of  measured  versus  predicted  data  for  approximately 
20  variables.  The  number  of  comparisons  for  each  variable  on  each  sam¬ 
pling  day  was  usually  between  7  and  10,  with  between  12  and  32  sampling 
dates  for  each  variable  for  each  year.  In  1981,  a  total  of  3,408  com¬ 
parisons  were  made;  in  1982,  the  number  was  3,251. 

29.  Because  of  the  large  amounts  of  measured  data  that  were  avail¬ 
able  for  model  confirmation,  adequate  judgment  of  the  total  model  was 
difficult  by  only  viewing  graphs.  A  change  in  a  particular  algorithm  or 
a  different  set  of  coefficients  may  have  improved  some  variables,  or 
improved  some  profiles  of  a  particular  variable,  while  making  others 
less  desirable.  For  those  cases,  statistical  analysis  for  all  compari¬ 
sons  was  helpful.  These  analyses  were  used  to  test  which  of  two  algo¬ 
rithms  for  a  particular  process  was  a  better  predictor  or  which  of  a 
number  of  sets  of  coefficients  produced  simulation  curves  that  most 
closely  corresponded  to  all  observed  data. 

30.  The  statistic  used  for  comparison  was  the  Reliability 
Index  (RI)  of  Leggett  and  Williams  (1981).  It  is  a  good  statistic  for 
aggregating  and  comparing  results  of  different  variables  because  it  does 
not  depend  on  whether  the  observed  or  predicted  value  is  greater;  also, 
it  is  scale  variant  against  additive  variation. 

31.  The  statistic  was  computed  for  each  variable  on  each  sampling 
day,  as  well  as  for  each  variable  over  depth  and  time.  An  average  value 
for  all  composite  RIs  was  calculated  to  give  the  "goodness"  for  the  en¬ 
tire  model.  Calculated  values  of  RI  could  range  from  1.0,  for  the  case 
of  perfect  prediction,  to  infinity.  If  all  comparisons  had  measured 
versus  predicted  values  within  a  factor  of  two  of  each  other,  the  value 
for  the  RI  would  be  2.0.  An  RI  of  10.0  signifies  that  the  differences 
between  measured  and  predicted  values  were  an  average  of  one  order  of 
magnitude  apart,  while  an  RI  of  100.0  signifies  that  the  comparisons 
were  two  orders  of  magnitude  apart. 

32.  Besides  comparing  the  concentration  of  predicted  versus  mea¬ 
sured  variables,  predicted  flux  values  were  compared  to  their  measured 
counterpart.  Data  were  available  on  algal  settling,  gross  production, 
and  oxygen  utilization  at  the  sediment-water  interface.  Settling  and 
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gross  production  values  were  measured  in  situ.  Oxygen  utilization 
values  were  measured  in  laboratory  studies  (Gunnison,  Chen,  and  Brannon 
1983).  Predicted  flux  values  were  not  checked  for  each  calibration 
simulation  because  the  flux  utility  is  used  after  a  CE-QUAL-R1  simula¬ 
tion  is  finished  (Wlosinski  1984). 
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PART  TV:  RESULTS  AND  DISCUSSION 
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Calibration  Simulations 

33.  Calibration  of  CE-QUAL-RI,  using  the  1981  Eau  Galle  data  set, 
required  a  number  of  phases  that  illustrated  the  difficulty  of  producing 
a  model  that  is  always  general  enough  to  simulate  all  reservoirs  while 
only  changing  model  coefficients.  Several  site-specific  modifications 
had  to  be  made  to  the  model  to  describe  adequately  some  of  the  charac¬ 
teristics  peculiar  to  Eau  Galle.  These  changes  are  discussed  in  the 
following  paragraphs. 

34.  Volume  of  the  Eau  Galle  Reservoir  was  calculated  from  sedi¬ 
ment  range  studies  performed  by  the  USAED,  St.  Paul,  in  1977.  The 
relationship  between  measured  volume  and  elevation  is  presented  in  Fig¬ 
ure  5  (curve  A).  Model  predicted  volumes  were  originally  computed  from 
the  area-elevation  equation. 
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Area  =  5 , 791 . EL 


1.954 


where 


EL  =  elevation  from  the  pool  bottom  (m) 

2 

Area  =  area  of  the  pool  at  a  particular  elevation  (m  ) 

The  coefficients  of  Equation  1  were  determined  from  the  sediment  studies 
using  regression  techniques.  Volume,  predicted  as  a  function  of  eleva¬ 
tion,  is  also  presented  in  Figure  5.  Because  the  predicted  volumes  were 
considered  unsatisfactory,  other  relationships  were  tried,  and  the  area- 
elevation  algorithm  was  replaced  by 

Area  =  4, 906. EL2  -  5, 707. EL  +  48,453  (2) 
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the  original  dimensions  describing  the  valley.  Figure  6  represents  con¬ 
tours  for  the  conservation  pool  which  were  obtained  from  the  1977  sedi¬ 
ment  range  studies.  The  area  below  279  m  represents,  approximately,  the 
borrow  area.  This  change  produced  satisfactory  temperature  predictions, 
having  an  average  RI  of  1.12. 


Figure  6.  Depth  contours  for  Eau  Galle  Reservoir 

37.  After  this  change  was  made,  the  CE-THERM-R1  data  set  was  then 
used  as  a  basis  for  the  CE-QUAL-R1  data  set.  The  first  set  of  predic¬ 
tions,  using  coefficients  from  the  literature  and  other  modeling  stud¬ 
ies,  was  considered  unsatisfactory.  The  average  RI  for  this  simulation 
was  6.5,  compared  to  2.8A  and  2.59  obtained  for  the  final  DeGray  Lake 
simulations  for  1979  and  1980,  respectively  (Wlosinski  and  Collins 
1985).  Subsequently,  a  number  of  simulations  were  made  during  which 
biological  and  chemical  coefficients  were  changed.  Although  predictions 
improved,  they  were  still  considered  unsatisfactory,  having  an  average 
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K1  of  5.03.  Evidence  from  these  simulations  pointed  toward  problems 
with  inflow,  outflow,  and  mixing  in  the  lower  layers  of  the  reservoir. 

38.  At  that  time,  conversations  with  personnel  of  the  I'SAED, 

St.  Paul,  led  to  the  discovery  that  the  old  river  and  creek  channels  may 
have  been  diverted  from  the  borrow  area.  Station  20,  where  the  vertical 
profiles  for  calibration  and  confirmation  data  were  measured,  was  in  the 
deepest  part  of  the  reservoir,  which  was  located  in  the  middle  of  the 
borrow  area.  Unlike  many  reservoirs  this  deep  area  was  not  next  to  the 
dam  (Figure  6),  in  the  old  river  channel.  In  a  typical  reservoir, 
cooler  water  entering  as  inflow  follows  the  old  thalweg  and  can  replace 
the  water  in  the  deepest  part  of  the  reservoir.  This  cooler  water  may 
also  be  lost  as  outflow  if  the  dam  is  equipped  for  low-level  withdrawal. 
Because  the  model  representation  befits  this  typical  reservoir,  changes 
in  pool  geometry  that  alter  the  normal  flow  pattern  would  necessitate 
model  changes  representing  flow. 

39.  Evidence  did  exist  to  show  that  construction  activities 
altered  these  normal  flow  patterns.  On  the  conservation  pool  borrow 
area  plan,  dated  20  May  1965  (USAED,  St.  Paul,  drawing  number 

M31a. 1  — L— 10/6) ,  a  note  was  made  that  the  Eau  Galle  River  and  Lohn  and 
I.ousv  Creeks  may  be  diverted  outside  the  borrow  area  limits.  This  is 
shown  in  Figure  7,  a  picture  taken  on  23  May  1968,  3  months  before  the 
dam  was  completed.  Point  A  marks  the  borrow  area;  B,  the  inlet  of  the 
Eau  dalle  River;  and  C,  the  channelized  outlet  of  Lohn  and  Lousy  Creeks. 
As  can  be  seen,  the  river  and  creeks  do  not  enter  the  borrow  area. 

Also,  a  mixing  study  performed  during  the  summer  of  1981  supported  this 
hypothesis.  During  that  study,  water  from  the  Eau  Galle  River  entered 
the  reservoir  as  an  underflow  at  station  40  and  continued  as  an  under¬ 
flow  to  station  30  (Figure  1).  It  then  spread  across  the  lake  as  an 
interflow  at  a  depth  from  3-4  m.*  It  was  therefore  concluded  that  the 
lower  elevations  of  the  borrow  area  were  isolated  from  inflow  and  out¬ 
flow  effects.  To  better  represent  this  isolation,  it  was  assumed  that 
the  lowest  5  m  of  the  pool  was  cut  off  to  inflowing  and  outflowing 


*  Personal  Communication,  February  1983,  Marc  Johnson,  US  Army  Engineer 
Waterways  Experiment  Station,  Vicksburg,  Miss. 
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Figure  7.  Eau  Caile  Reservoir,  3  months  before  dam  completion 
(A  is  located  in  the  borrow  pit,  B  is  located  at  the  inflow  of 
the  Eau  Galle  River,  and  C  is  located  at  the  inflow  of  Lohn 

and  Lousy  Creeks) 

waters.  This  step  aided  model  predictions,  but  problem  areas  still 
remained . 

40.  It  appeared  at  this  time  that  some  of  the  variable  concentra¬ 
tions  in  the  metalimnion,  during  the  summer  months,  were  reflecting  too 
much  input  from  inflowing  waters,  whereas  the  predictions  in  the  epilim- 
nion  would  be  better  if  some  of  the  inflowing  water  was  mixed  with  sur¬ 
face  waters.  In  the  original  model,  inflow  placement  was  dependent  only 
on  the  density  of  inflowing  water  compared  to  the  density  within  the 
pool.  There  is  no  mixing  or  entrainment  with  other  layers  in  the  reser¬ 
voir  unless  density  is  nearly  the  same.  As  can  be  seen  in  Figure  1, 
reservoir  tributaries,  especially  Lohn  and  Lousy  Creeks,  are  relatively 

3 

small,  and  since  flows  are  usually  less  than  1  m  /sec,  some  mixing  is 
quite  reasonable.  To  test  this  hypothesis,  an  algorithm  that  mixed 


inflowing  water  with  the  surface  layer  as  a  function  of  wind  was  added 
to  the  model.  The  algorithm  is  depicted  in  Figure  8.  The  fraction  of 
inflowing  water  that  mixes  with  the  surface  increases  with  an  increase 
in  wind  or  a  decrease  in  inflow  discharge.  This  model  change  was  tempo¬ 
rary,  being  made  only  for  the  Eau  Galle  simulations.  Although  improve¬ 
ments  were  made  concerning  predictions,  results  were  still  unsatisfac¬ 
tory.  Problems  at  this  time,  though,  did  not  appear  to  be  due  to  mixing 
or  flow  problems,  but  to  problems  with  the  biological  portion  of  the 
model . 


WIND,  km/hr 

Figure  8.  The  relationship  between  the  amount  of 
wind  and  the  amount  of  inflowing  water  mixing  in 
the  surface  layer 

41.  Orthophosphate  phosphorus  was  measured  on  7  April  1981,  with 
values  ranging  in  the  water  column  from  0.008  to  0.011  mg/Z.  It  was  not 
possible  to  predict  the  diatom  bloom  of  169  ug/i-  chlorophyll-a  measured 
on  21  April  when  using  the  measured  orthophosphate  values.  Even  if  all 
of  the  orthophosphate  measured  was  converted  to  algae,  with  no  algal 


losses  to  sinking,  respiration,  zooplankton  ingestion,  or  other  factors, 
the  level  of  algae  would  not  be  as  high  as  measured  on  21  April.  During 
this  2-week  period,  flow  levels  remained  low,  as  did  measured  concentra¬ 
tions  of  orthophosphate  in  the  inflow.  On  4  April,  3  days  before  the 

start  of  simulation,  1.46  in.  of  rainfall  caused  the  flow  of  the  Eau 

3  3 

Galle  River  to  rise  from  a  normal  0.2-0. 9  m  /sec  to  6.7  m  /sec.  This 
rise,  along  with  a  probable  increase  in  orthophosphate  concentrations 
(Ashby  and  James  1985) ,  may  have  supplied  the  phosphorus  needed  for 
growth.  The  diatoms  may  have  increased  their  internal  phosphorus  con¬ 
centrations  by  way  of  luxury  uptake  (Rhee  1974,  Collins  1980).  Luxury 
uptake  is  not  modeled  in  CE-QUAL-R1,  and  to  do  so  would  make  the  model 
more  complex  than  is  needed  for  the  Corps.  To  overcome  this  problem, 
the  initial  conditions  for  orthophosphate  were  increased  to  a  value  be¬ 
tween  measured  orthophosphate  and  total  phosphorus.  For  the  1982  con¬ 
firmation  simulations,  the  measured  values  (near  0.08  mg/£)  for  ortho¬ 
phosphate  were  used. 

42.  To  help  locate  other  problem  areas,  equivalent  changes  were 
made  to  the  flux  model,  with  flux  predictions  then  being  checked.  Pre¬ 
dictions  for  which  data  were  available  (algal  s  tting,  gross  production, 
and  oxygen  utilization  at  the  sediment-water  interlace)  appeared  reason¬ 
able.  Although  most  of  the  biological  variables  are  interconnected, 
with  changes  to  one  variable  affecting  the  predicted  results  of  other 
variables,  it  appeared  from  studying  the  results  that  most  of  the  pro¬ 
blems  were  due  to  an  excess  of  nitrogen  in  the  system  when  anaerobic 
conditions  existed,  and  to  the  zooplankton  compartment  not  responding 
properly  to  changes  in  food  concentrations.  Denitrification  can  occur 
intensely  in  anaerobic  environments,  a  process  by  which  nitrite  is 
reduced  to  elemental  nitrogen,  which  can  then  be  lost  from  the  system. 
This  process  was  permanently  added  to  the  model  after  simulations  showed 
improved  predictions. 

43.  To  improve  the  zooplankton-algae  relationship,  a  threshold 
value  at  which  grazing  commences  was  added  to  the  zooplankton  ingestion 
algorithm.  Experimental  evidence  on  a  grazing  threshold  was  supplied  by 
Parsons,  LaBrasseur,  and  Fulton  (1967)  and  McAllister  (1970).  However, 
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in  a  review  of  simulation  modeling  of  zooplankton,  Leidy  and  Ploskey 
(1980)  mention  that  the  need  for  a  grazing  threshold  may  stem  from  inap¬ 
propriate  assumptions  or  our  ignorance  of  the  grazing  dynamics  of  zoo¬ 
plankton.  This  may  very  well  be  true,  but  in  the  case  of  CE-QUAL-R1 , 
zooplankton  are  modeled  mainly  because  of  their  impact  on  other  vari¬ 
ables  that  are  of  greater  environmental  impact,  namely  oxygen  and  algae. 
Therefore,  to  restrain  additional  complexity,  only  the  threshold  for 
grazing  was  added. 

44.  Other  changes  to  the  model  that  were  made  at  this  time  in¬ 
cluded  a  minor  correction  to  the  adsorption  algorithm  and  simplifica¬ 
tions  to  the  anaerobic  portion  of  the  model.  As  with  many  other  chemi¬ 
cal  and  biological  rates  in  the  model,  many  of  the  anaerobic  processes 
were  temperature  dependent.  Unfortunately,  very  little  information  is 
available  in  the  literature  concerning  temperature  dependence.  Because 
the  temperature  coefficients  used  for  both  the  DeGray  Lake  and  Eau  Galle 
Reservoir  simulations  were  such  that  the  temperature  range  during  the 
year  would  have  no  effect  on  the  rates  of  anaerobic  processes,  it  was 
recommended  that  the  temperature  coefficients  be  removed.  This  change 
reduces  the  number  of  coefficients  by  78,  and  thus  saves  time  in  data 
set  preparation. 

45.  After  these  changes  were  made  and  a  number  of  calibration 
simulations  were  made  to  fine-tune  the  model,  comparisons  of  predicted 
to  measured  values  were  much  more  acceptable,  giving  an  overall  RI  of 
2.63.  A  confirmation  simulation,  using  the  same  coefficients  in  the 
1982  data  set  as  used  for  1981,  also  produced  acceptable  results,  giving 
an  overall  RI  of  2.71.  However,  the  fluxes  for  the  confirmation  were 
considered  unsatisfactory.  Predictions  of  gross  production  were  nearly 
an  order  of  magnitude  lower  than  measured  values,  algal  settling  predic¬ 
tions  were  nearly  an  order  of  magnitude  low,  and  sediment  oxygen  utili¬ 
zation  was  nearly  an  order  of  magnitude  high,  consuming  over  80  percent 
of  all  oxygen  used  for  respiration.  In  addition,  algal  respiration  was 
only  8  percent  of  gross  production,  an  extremely  low  value.  Thus,  it 
appeared  that  acceptable  predictions  were  being  made  for  the  wrong  rea¬ 
sons.  After  another  set  of  calibration  simulations  were  made,  while 
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adjusting  coefficients  dealing  with  the  above  rates  and  checking  flux 
predictions,  acceptable  results  were  again  obtained.  The  final  overall 
RI  for  the  1981  data  set  was  2.57,  and  for  the  1982  confirmation  data 
set,  2.62.  A  discussion  of  the  flux  predictions  for  the  final  simula¬ 
tions  is  provided  in  the  last  section  of  Part  IV. 

46.  Statistical  results  for  each  variable  from  the  final  calibra¬ 
tion  simulation  are  presented  in  Table  1  (see  Wlosinski  1984  for  a 
discussion  of  statistics).  Figure  9  presents  four  graphs  (representing 
different  dates)  of  each  of  19  variables,  comparing  predictions  (solid 
lines)  to  measured  values  (dots).  The  graph  for  the  date  with  the  poor¬ 
est  RI  is  included  for  each  variable.  Values  for  the  Rl  ranged  from 
1.07  to  4.85.  In  general,  variables  whose  concentration  normally  ranges 
over  more  than  one  order  of  magnitude  had  higher  RI  values.  For  exam¬ 
ple,  orthophosphate  phosphorus  values  ranged  from  0.001  to  over  1.0  mg/f. 
and  had  an  RI  of  4.84,  whereas  total  dissolved  solids,  which  ranged  from 
168  to  225  mg/) l,  had  an  RI  of  1.24.  Predictions  of  temperature,  which 
are  extremely  important  because  of  the  effect  of  temperature  on  most 
biological  processes,  were  considered  excellent,  having  an  RI  of  1.09. 
The  timing  of  initial  stratification  and  of  fall  overturn  was  good,  as 
was  the  shape  of  the  temperature  profile. 

47.  Most  of  the  major  dynamics  of  other  variables  were  predicted 
correctly.  Chlorophyll  predictions  in  the  epilimnion  were  usually 
within  a  factor  of  two.  Hypolimnetic  predictions  were  not  as  good,  but 
were  still  considered  acceptable.  Even  though  the  graphics  and  statis¬ 
tics  were  performed  on  the  summation  of  the  predicted  values  for  the 
three  algal  compartments,  individual  comparisons  in  the  epilimnion  can 
be  made  to  assess  the  model's  ability  to  predict  seasonal  succession  of 
major  groups.  As  stated  earlier,  Barko  et  al.  (1985)  found  that  algae 
were  dominated  by  the  Bacillariophyceae  (diatoms),  Cyanophyta  (blue- 
greens),  and  Pyrrophyta  (dinof lagellates) .  Using  data  supplied  in  their 
paper,  graphs  representing  seasonal  succession  of  the  three  major  com¬ 
partments  were  made  for  the  epilimnion  (Figure  10) .  The  solid  line 
represents  measured  data,  and  the  dots  represent  model  predictions.  It 
is  apparent  that  the  model  reasonably  predicted  the  seasonal  succession 
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Figure  9.  Comparison  of  predicted  values  to  measured  data  for  19  variables — final  Eau  Galle 
calibration  simulation,  1981.  The  solid  line  represents  model  predictions;  dots  represent 

measured  values  (Sheet  1  of  10) 
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of  the  dominant  algal  groups  measured  in  Eau  Galle  Reservoir  in  1981. 

48.  Most  of  the  predicted  dissolved  silica  values  were  above  mea¬ 
sured  values,  averaging  9.19  versus  6.05  mg/ l,  respectively.  The  prob¬ 
lem  may  be  due  parti}'  to  the  model's  constant  stoichiometric  equivalent 
between  algae  and  silica,  which  can  vary  in  the  field,  and  partly  to 
silica  adsorption  and  sedimentation  in  reservoirs,  a  process  not  in¬ 
cluded  in  the  model.  Silica  adsorption  should  be  considered  for  future 
model  applications  if  diatom  dynamics  are  important  for  a  particular 
study.  The  major  difference  of  both  ammonia-nitrogen  and  phosphate- 
phosphorus  concentrations  was  spatial,  in  that  concentration  differences 
of  over  two  orders  of  magnitude  between  the  surface  and  bottom  waters 
occurred.  Acceptable  predictions  of  these  differences  were  made. 

49.  The  dynamics  of  nitrite  plus  nitrate  were  more  difficult  to 
predict,  even  with  the  process  of  denitrification  added  to  the  model. 
Early  in  the  calibration  process,  problems  occurred  in  predicting  the 
maxima  that  often  occurred  at  3-6  m  (e.g.  Figure  9,  sheet  5,  date 
810630),  a  problem  that  was  later  solved  by  not  allowing  inflow  into  the 
layers  representing  the  borrow  area.  However,  on  14  July  the  bulge  dis¬ 
appeared,  only  to  recur  on  28  July.  This  disappearance  was  not  predicted 
by  the  model.  The  low  values  in  the  epilimnion  and  hypolimnion  during 
stratified  periods  were  predicted  correctly,  but  for  different  reasons. 

In  the  epilimnion,  the  loss  from  the  compartment  was  due  to  the  utiliza¬ 
tion  of  the  nutrient  by  way  of  photosynthesis,  whereas  the  loss  in  the 
hypolimnion  was  due  to  denitrification.  On  14  July,  virtually  no  oxygen 
was  measured  at  station  20  below  3  m,  whereas  the  model  predicted  up  to 

3  mg/ l  between  3  and  5  m.  Since  denitrification  occurs  only  under 
anaerobic  conditions,  the  positive  oxygen  predictions  between  3  and  5  m 
probably  caused  the  problem  with  the  nitrite-nitrate  compartment. 

50.  Although  most  of  the  oxygen  predictions  were  quite  close  to 
measured  values,  the  scenario  above  illustrates  the  problems  that  can 
occur  if  the  predictions  are  not  exact.  In  this  case  the  problem  was 
due,  to  a  large  measure,  to  the  one-dimensional  assumption  and  the  mix¬ 
ing  of  inflowing  waters.  At  3  m  below  the  surface  at  station  20,  where 
the  confirmation  data  were  measured,  the  concentration  of  dissolved 
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oxygen  was  0.5  mg/i  and  the  concentration  of  nitrite-nitrate  was 
0.003  mg /i.  At  station  30  (Figure  1),  at  the  same  depth  and  time,  oxy¬ 
gen  was  measured  as  6.8  mg/£.  and  nitrite-nitrate  as  0.296  mg/i.  Thus, 
denitrification  could  have  occurred  at  station  20,  but  not  at  30.  This 
was  probably  due  to  the  influence  of  inflowing  waters,  which  affect  sta¬ 
tion  30  more  than  station  20.*  On  July  14,  oxygen  concentrations  in  the 
inflowing  waters  of  the  Eau  Galle  River  and  French  Creek  were  near 
6.4  mg/£,  and  nitrite-nitrate  concentrations  were  near  0.8  mg/i.  Tem¬ 
peratures  for  the  two  tributaries  were  18.0  and  14.5°  C,  respectively, 
whereas  the  temperature  at  station  30  ranged  from  24.7  to  21.7°  C.  In¬ 
flowing  water  probably  entered  the  reservoir  as  an  underflow  and  fol¬ 
lowed  the  old  thalweg  to  station  30. 

51.  Because  nitrite-nitrate  values  at  station  20  ranged  from 
undetectable  limits  to  0.03  mg/i,  with  values  below  the  dam  at  0.5  mg/i, 
it  is  suggested  that  under  some  conditions  some  of  the  inflowing  water 
may  have  followed  the  thalweg  to,  and  through,  the  outflow  port.  For 
these  same  reasons,  care  must  be  taken  when  using  outflow  concentrations 
predicted  by  the  model  if  the  deepest  part  of  the  reservoir  is  not  near 
the  outlet  ports.  Conditions  in  the  littoral  zone  may  then  dominate  the 
outflow. 

52.  Two  other  variables  exhibited  concentration  differences  be¬ 
tween  the  midlayers  and  the  epilimnion  and  hypolimnion,  this  time  in  the 
form  of  minima.  Both  total  organic  carbon  and  dissolved  organic  carbon 
had  lower  concentrations  between  3  and  6  m  than  found  in  the  epilimnion 
or  hypolimnion  and,  again,  lower  values  in  the  inflow  that  were  placed 
in  the  metalimnion  could  help  explain  the  dynamics.  The  model  predicted 
these  minima. 

53.  There  was  one  other  event  that  was  not  predicted  at  all  by 
the  model.  Sulfate,  the  concentrations  of  which  were  normal  through 
August  for  a  system  with  a  clinograde  oxygen  curve,  fell  to  near  zero 
for  four  sampling  periods  in  September  and  October.  Unlike  the  problem 


*  Personal  Communication,  November  1983,  J.  H.  Carroll,  US  Army 
Engineer  Waterways  Experiment  Station,  Vicksburg,  Miss. 
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with  nitrite  plus  nitrate,  these  low  values  were  found  at  all  pool  sta¬ 
tions,  including  station  40,  as  well  as  the  downstream  station.  Even 
though  some  of  the  inflowing  concentrations  for  this  period  were  low, 
others  were  as  high  as  14.5  mg/{,.  Although  field  personnel  noticed  the 
low  concentrations  and  authenticated  the  results,*  they  were  unable  to 
establish  a  cause  for  the  phenomenon.  Sulfate  dynamics  during  other 
months,  as  well  as  the  dynamics  of  other  anaerobic  variables,  were 
predicted  satisfactorily. 

54.  Macrophyte  biomass  measurements  were  taken  in  1981  (Filbin 
and  Barko  1985).  Because  their  measurements  were  taken  on  an  aerial 
basis,  the  RI  statistic  was  not  used  for  comparing  measured  and  pre¬ 
dicted  values.  Filbin  and  Barko  reported  that  the  macrophyte  biomass 
and  associated  epiphytes  occupied  about  17  percent  of  the  surface  area 
of  Eau  Galle.  This  figure  was  used  to  convert  the  predicted  mass  from 
units  of  grams  per  reservoir  to  the  measured  units  of  grams  per  square 
metre.  A  comparison  of  predicted  versus  measured  values  is  presented  in 
Figure  11.  Considering  that  the  model  lumps  all  macrophyte  and  epiphyte 
species  together,  predictions  are  considered  acceptable. 
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Confirmation  Simulation 


55.  Statistical  results  from  the  final  confirmation  simulation  of 
1982  are  presented  in  Table  2.  Comparisons  of  measured  versus  predicted 
values  for  selected  dates  are  presented  in  Figure  12.  As  for  1981,  the 
graph  representing  the  date  with  the  poorest  RI  is  included  for  each 
variable.  The  average  RI  for  all  variables  was  a  satisfactory  2.62, 
with  a  range  from  1.06  to  4.85.  The  RI  for  temperature  was  1.14, 
although  on  a  few  dates  the  predicted  mixed  layer  was  approximately  1  m 
too  deep.  This  was  probably  due  to  the  wind  difference  between  the  site 
and  Eau  Claire  (see  paragraph  20),  since  there  was  an  appreciable  differ¬ 
ence  in  temperature  predictions  depending  on  the  station  where  meteoro¬ 
logical  data  were  measured.  The  model  failed  to  accurately  predict  an 


*  Personal  Communication,  J.  H.  Carroll,  op.  cit. 
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1981 

Figure  11.  Measured  values  (circles)  and  predicted 
values  (triangles)  for  macrophyte  biomass  for  Eau 
Galle  Reservoir  in  1981 

algal  bloom  that  was  measured  on  18  May,  which  probably  also  accounted 
for  the  high  dissolved  silica  prediction  in  the  epilimnion  when  the  mea¬ 
sured  values  were  quite  low.  Most  of  the  other  algal  predictions  were 
reasonable,  with  a  predicted  average  of  19.4  Ug/i  versus  a  measured 
average  of  23.5  Ug/l.  As  in  1981,  most  silica  predictions  were  slightly 
high  and  were  probably  due  to  the  lack  of  silica  dynamics  in  the  model. 

56.  A  slight  error  in  predicting  temperature  can  magnify  the 
errors  in  predicting  other  variables.  This  occurred  on  5  October,  when 
the  measured  temperature  profile  indicated  that  fall  turnover  had  not 
yet  occurred,  and  anaerobic  conditions  existed  8  m  below  the  surface. 
Complete  mixing  was  predicted  before  this  date,  with  the  effect  that 
concentrations  of  ammonia,  nitrite  plus  nitrate  nitrogen,  suspended 
solids,  total  manganese,  total  iron,  sulfate,  reduced  manganese,  and 
reduced  iron  were  poorly  predicted  for  that  date.  The  dynamics  of  these 
variables  were  predicted  correctly;  it  was  the  timing  that  was  incor¬ 
rect.  In  actual  applications  of  the  model,  it  is  much  more  important  to 
predict  the  dynamics  of  the  system  rather  than  the  timing.  This  is 
because  the  timing  is  due,  to  a  great  extent,  to  meteorological 
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conditions,  which  themselves  cannot  be  accurately  transferred  from  the 
meteorological  station  to  the  site. 

57.  The  maxima  for  nitrite  plus  nitrate  and  the  minima  for  dis¬ 
solved  organic  carbon  and  total  organic  carbon  occurred  in  1982  similar 
to  1981.  Although  these  phenomena  were  again  predicted  by  the  model, 
the  concentrations  were  not  as  accurate  as  the  1981  predictions.  Unlike 
1981,  sulfate  concentrations  in  the  epilimnion  remained  around  10  mg/i 
through  the  summer  and  fall,  with  excellent  predictions  being  made  by 
the  model.  For  sulfate,  the  RI  for  1982  was  1.86,  compared  to  3.41  for 
1981. 

58.  The  above  results  have  a  bearing  on  possible  uses  of  the 
model.  in  the  case  of  pre impoundment  studies,  because  unsatisfactory 
results  were  obtained  from  the  initial  calibration  simulation,  results 
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would  be  questionable  unless  real  data  were  available  from  similar 
impoundments  nearby.  Johnson  and  Ford  (1981)  confirmed  the  model 
(CE-THERM-R1 )  in  this  manner  on  DeGray  Lake  and  Lake  Greeson  in  Arkan¬ 
sas,  and  obtained  excellent  results.  To  date,  the  same  type  of  con¬ 
firmation  exercise  has  not  been  made  for  the  full  CE-QUAL-R1  model.  In 
addition,  the  range  of  results  obtained  when  using  wind  data  from  two 
different  stations  for  the  Eau  Galle  simulation  leads  to  the  recommenda¬ 
tion  that  both  pre impoundment  and  post impoundment  studies  include  simu¬ 
lations  using  meteorological  data  for  a  number  of  different  years.  Fur¬ 
thermore,  because  of  differences  between  reservoirs,  minor  changes  to 
the  computer  code  may  have  to  be  made  to  address  a  specific  problem. 


Time  Steps 


59.  One  of  the  features  of  CE-QUAL-R1  is  a  variable  time  step, 
allowing  the  user  to  choose  from  1-,  2-,  3-,  4-,  6-,  8-,  12-  or  24-hr 
periods.  Shorter  time  steps  should  show  some  of  the  dynamics  occurring 
during  a  1-day  time  period,  although  simulations  would  be  more  expen¬ 
sive.  All  of  the  simulations  for  evaluating  the  model  while  using  data 
collected  on  DeGray  Lake,  as  well  as  other  work  reported  here,  used 
24-hr  time  steps.  To  test  the  model  using  periods  other  than  24  hr,  a 


number  of  simulations  were  made  using  a  3-hr  time  step.  The  first  simu¬ 
lation  used  the  1981  calibration  data  set  from  Eau  Galle,  without  any 
changes  to  coefficients.  In  addition,  the  same  set  of  driving  variables 
was  used.  The  model  simply  interpolated  the  24-hr  data  to  supply  values 
for  3-hr  periods.  It  would  have  been  more  appropriate  to  supply  data 
measured  at  3-hr  intervals. 

60.  The  results  were  not  at  all  satisfactory,  indicating  that  the 
same  set  of  coefficients  cannot  be  used  with  different  time  steps.  This 
problem  is  probably  due  to  the  nonlinearity  of  some  equations.  After  a 
number  of  calibration  simulations  were  made,  reasonable  results  were  ob¬ 
tained.  To  produce  satisfactory  results  for  the  calibration  simulation 
using  3-hr  time  steps,  eight  coefficients  had  to  be  changed.  These  were 
the  mixing  coefficient  due  to  wind,  the  light  saturation  and  maximum 
gross  production  coefficients  for  the  three  algal  compartments,  and  the 
zooplankton  ingestion  coefficient.  Although  the  predictions  for  macro¬ 
phytes  were  not  checked,  it  seems  reasonable  that  for  different  time 
steps,  the  macrophyte  gross  production  and  light  saturation  coefficients 
should  also  be  changed.  The  RI  for  the  final  3-hr  simulation  was  2.69, 
although  only  a  few  calibration  runs  were  made  compared  to  the  24-hr 
time  step.  It  is  recommended  that  once  the  model  has  been  calibrated 
for  a  set  of  coefficients  and  time  step,  the  time  step  should  not  be 
changed. 

Flux  Predictions 

61.  Wlosinski  (1979),  Scavia  (1980),  Chapra  et  al.  (1983),  and 
Collins  and  Wlosinski  (1984)  have  shown  the  need  to  compare  predicted 
and  measured  flux  values  in  addition  to  mass  or  concentration  values. 
Evaluating  models  by  concentrating  solely  on  the  comparison  of  measured 
versus  predicted  concentrations  can  produce  a  model  that  predicts  rea¬ 
sonable  concentrations  for  the  wrong  reasons.  As  part  of  the  Eau  Galle 
Reservoir  study,  a  number  of  process  rates  were  measured  which  could  be 
compared  to  model  predictions.  These  included  algal  productivity,  sedi¬ 
ment  oxygen  demand,  and  sedimentation  rates  of  algae. 


PRODUCTION 


62.  Productivity  rates  vary  extensively  over  a  24-hr  period  and 
from  day  to  day,  and  are  dependent  on  local  conditions  such  as  light, 
temperature,  and  the  concentrations  of  nutrients  and  algae.  Neverthe¬ 
less,  comparisons  should  show  rates  that  are  similar.  Production  rates 
were  measured  from  noon  until  2  pm  at  discrete  intervals  between  the 
surface  and  3  m  on  a  biweekly  sampling  schedule.  The  measured  values  of 
mg  02/m2/hr  were  converted  to  g  O^/mVday,  assuming  a  10-hr  photosynthe¬ 
tic  day,  for  comparison  with  predicted  values.  The  comparison,  shown  in 
Figure  13,  used  predictions  from  an  approximate  depth  of  2  m.  It  must 
be  remembered  that  algal  processes  and  rate  values  used  in  the  model  are 
based  upon  average  daily  values.  Studies  have  revealed  that  photosyn¬ 
thetic  capacity,  cell  division,  nutrient  uptake,  respiration,  and  graz¬ 
ing  vary  according  to  algal  circadian  rhythms  (Prezelin  et  al.  1977; 


Chisholm,  Azam,  and  Eppley  1978;  and  Chisholm  1981).  Considering  the 
variability  in  measured  data  at  different  stations  and  the  fact  that  the 
measured  data  are  for  a  much  shorter  period  than  the  24-hr  time  step 
used,  predictions  are  quite  reasonable. 

63.  Sedimentation  rates  of  algae  were  also  available.  Sediment 


Figure  13.  Algal  net  production  rates  for  Eau  Galle  Reservoir,  1981 
station  20.  Circles  represent  observed  rates  from  0-3  m;  the  solid 
line  represents  predicted  rates  at  approximately  2  m 


traps  were  deployed  at  4  and  8.5  m  below  the  water  surface  and  were  re¬ 
trieved  after  approximately  2  weeks.  Direct  comparison  of  predicted  and 
measured  values  was  not  possible  because  of  the  variable  layer  scheme 
employed  in  the  model.  The  model  predicts  values  on  a  per-layer  basis, 
and  because  the  thickness  of  each  layer  is  variable,  the  4-  or  8.5-m 
depth  contour  may  be  located  in  different  layers  during  different  time 
steps.  For  comparative  purposes,  the  predictions  from  layer  4  (Fig¬ 
ure  14)  and  layer  8  (Figure  15)  were  used.  In  general,  the  predictions 
were  more  uniform  than  the  measured  data,  and  except  in  June  at  4  m,  the 
model  predictions  were  lower  than  measured  values.  Part  of  the  differ¬ 
ence  may  be  due  to  the  fact  that  as  algae  settle  in  the  model,  their 
mass  is  lost  to  the  nutrients,  labile  dissolved  organic  matter,  and 
detritus  through  respiration  and  mortality,  whereas  chlorophyll  a  in  the 
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system  decomposes  more  slowly  and  may  be  measured  at  deeper  depths  even 
though  the  cells  are  not  viable.  Most  of  the  values  are  close  enough  to 
be  considered  reasonable,  especially  when  compared  to  earlier  predic¬ 
tions  that  were  an  order  of  magnitude  different  than  measured  values. 

64.  Sediment  oxygen  demand  was  measured  in  the  laboratory  on 
sediment  collected  from  Eau  Galle  Reservoir  (Gunnison,  Chen,  and  Brannon 
1983).  The  dissolved  oxygen  depletion  rate  was  176  mg/m2/day  at  20°  C, 
which  was  converted  to  modeled  units  of  32  kg/layer /day .  In  making  this 
conversion,  predictions  from  layer  4  were  used,  which  has  an  area  of 
183,000  m2.  Predicted  sediment  oxygen  demand  varies  extensively  over 
the  year,  but  usually  ranges  from  near  0  to  120  kg/layer/day  (Fig¬ 


ure  16).  The  major  differences  seen  are  a  function  of  temperature. 

65.  A  detailed  evaluation  of  all  processes  and  constituents  in¬ 


volved  in  dissolved  oxygen  production  and  consumption  indicates  that 
algal  production  (Figure  17)  and  macrophyte  production  (Figure  18) 
evolved  more  oxygen  which  was  consumed  by  chemical  and  biological  pro¬ 
cesses.  The  consuming  processes  include:  algae,  macrophyte,  zooplank¬ 
ton,  and  fish  respiration;  decomposition  of  sediment,  detritus,  and 
dissolved  organic  matter;  ammonia  decay;  and  anaerobic  oxidation 
(Table  3).  Correct  predictions  of  anaerobic  conditions  were  made  de¬ 
spite  the  prediction  of  greater  oxygen  production  versus  utilization. 


Figure  14.  Algal  dry  weight  sedimentation  rates  for  Eau 
Galle  Reservoir,  1981,  station  20.  Circles  represent 
observations  at  4  m,  and  squares  represent  predicted 

rates  for  layer  4 
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Figure  19.  Predicted  exchange  of  oxygen  at  the  air-water 
interface  for  Eau  Galle  Reservoir,  1981.  Negative  fluxes 
represent  oxygen  leaving  the  system 


PART  V:  SUMMARY  AND  CONCLUSIONS 


66.  Simulation  of  Eau  Galle  Reservoir  was  the  final  step,  during 
the  EWQOS  Program,  in  developing  the  water  quality  model  CE-QUAL-R1. 
Evaluation  has  included  tests  of  the  code  to  ensure  that  computer  pro¬ 
gramming  was  correct,  as  well  as  comparisons  of  model  predictions  to 
field  measured  values.  Tests  of  the  code  included  evaluations  of  the 
stability  of  predictions,  conservation  of  mass,  time  step  comparisons, 
entries  of  initial  values,  comparisons  using  different  driving  vari¬ 
ables,  and  a  check  of  equation  dimensionality  (Wlosinski  and  Collins 
1985)  . 

67.  Two  data  sets  (DeGray  Lake,  Ark.,  and  Eau  Galle  Reservoir, 
Wis.)  were  collected  specifically  to  evaluate  model  predictions.  The 
two  reservoirs  differed  markedly  in  size,  shape,  location,  withdrawal 
structure,  and  biological  attributes.  Evaluation  of  model  predictions 
and  field  measured  values  included  both  graphical  and  statistical  com¬ 
parisons.  The  confirmation  simulation  for  DeGray  Lake  had  an  average  RI 
of  2.59  (Wlosinski  and  Collins  1985),  and  for  Eau  Galle  Reservoir,  2.62. 
This  is  a  major  step  over  the  once  "accepted"  modeling  phrase  of  "accep¬ 
table  within  an  order  of  magnitude,"  which  was  used  for  algae  and  nu¬ 
trients.  All  variables  for  the  confirmation  exercises  had  RI  values 
under  5.0.  In  addition  to  comparisons  of  the  concentrations  of  vari¬ 
ables,  a  number  of  flux  values  were  also  investigated  for  the  two  reser¬ 
voirs,  and  all  were  considered  acceptable.  This  helps  ensure  that  var¬ 
iable  concentrations  were  predicted  for  the  correct  reasons.  This  is 
necessary  when  the  model  is  to  be  used  for  comparing  different  engineer¬ 
ing  strategies. 
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Results  from  the  1981  Final  Calibration  Simulation  of  Eau  Galle  Reservoir 


Total 


Calibration  Simulation  of  Eau  Galle  Reservoir 


Table  3 

Eau  Galle  1981  Oxygen  Budget* 


_ Process _ 

Algal  photosynthesis 
Macrophyte  photosynthesis 
Algal  respiration 
Macrophyte  respiration 
Zooplankton  respiration 
Fish  respiration 
Ammonia  decay 
Detritus  decomposition 
Sediment  decomposition 
Labile  DOM  decomposition 
Refractory  DOM  decomposition 
Anaerobic  oxidation 


Oxygen  Production  Oxygen  Utilizat 
0.1945  x  106 
0.2135  x  106 

0.4734  x  105 
0.8863  x  io5 
0.2815  x  105 
0.6555  x  10A 
0.4152  x  ioA 
0.5678  x  io4 
0.1218  x  io6 
0.2373  x  io5 
0.6165  x  io4 
0.3762  x  io3 


Total  oxygen  produced  0.4080  x  10^ 

Total  oxygen  utilized 


0.3325  x  106 
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ID 

EAU 

GALLE 

1981  JUN 

4  84 

#1 

FILES 

24 

353 

M2 

STPAL 

81 

113 

.2 

-6.1 

-9.8 

980.5 

11.8 

M2 

STPAL 

81 

114 

1  .  0 

-4.2 

-8.1 

982.5 

15.5 

U'2 

STPAL 

81 

115 

.4 

-12.2 

-18.8 

989.9 

21.5 

M2 

STPAL 

81 

116 

.1 

-13.3 

-19.2 

999.2 

14.6 

M2 

STPAL 

81 

117 

.5 

-5.3 

-12.6 

995.5 

19.4 

M2 

STPAL 

81 

118 

.2 

-.8 

-7.8 

991.8 

10.9 

M2 

STPAL 

81 

119 

.0 

-1.4 

-6.7 

982.1 

8.3 

M2 

STPAL 

81 

120 

.1 

-3.3 

-7.6 

987.1 

9.7 

M2 

STPAL 

81 

121 

.6 

-1.3 

-6.7 

990.6 

12.0 

M2 

STPAL 

81 

122 

.2 

.1 

-6.6 

987.7 

9.7 

M2 

STPAL 

81 

123 

.6 

.2 

-6 . 9 

983.5 

7.9 

M2 

STPAL 

81 

124 

.2 

3.5 

-3.9 

977.2 

11.8 

M2 

STPAL 

81 

125 

.6 

4.5 

-1.9 

967.7 

15.7 

W2 

STPAL 

81 

126 

1.0 

-1.3 

-6 . 9 

971.8 

23.6 

M2 

STPAL 

81 

127 

.7 

-7.3 

-12.5 

978.4 

16.2 

M2 

STPAL 

81 

128 

.5 

-10.2 

-17.2 

983.1 

18.8 

M2 

STPAL 

81 

129 

.1 

-10.9 

-17.6 

997.2 

13.0 

M2 

STPAL 

81 

130 

.2 

-9.2 

-16 . 0 

10  00.8 

11.6 

_M2 

STPAL 

81 

131 

1 .  0 

-4.0 

-11.8 

985.4 

18 . 3 

M2 

STPAL 

81 

2  1 

.7 

-8.3 

-14.0 

971.8 

28.7 

M2 

STPAL 

81 

2  2 

.3 

-19.4 

-26 . 1 

981.6 

23.4 

M2 

STPAL 

81 

2  3 

.4 

-16.8 

-23.2 

985.7 

18.8 

A4 


142 

STPAL 

al 

2  4 

.  1 

-16.5 

-22.8 

991.7 

18.8 

142 

STPAL 

81 

2  5 

.  9 

-8.3 

-13.6 

982 . 2 

18.8 

142 

STPAL 

81 

2  6 

.7 

-5.6 

-10.5 

982 . 9 

15.7 

142 

STPAL 

81 

2  7 

1  .  0 

-6.7 

-10.5 

971.6 

22.7 

142 

STPAL 

81 

2  8 

.7 

-12.4 

-18.0 

978.6 

14.1 

142 

STPAL 

81 

2  9 

.7 

-16 . 1 

-21.7 

980 . 3 

13.7 

142 

STPAL 

81 

210 

1.0 

-18 . 9 

-25.2 

983. 1 

24.5 

142 

STPAL 

81 

211 

.  1 

-23.3 

-29.7 

987 . 0 

16.2 

142 

STPAL 

81 

212 

.  5 

-15.1 

-20.9 

996.2 

10.2 

142 

STPAL 

81 

213 

.5 

-6.3 

-11.5 

999.2 

13.4 

142 

STPAL 

81 

214 

.8 

-1  .  9 

-6.2 

996.4 

15.7 

142 

STPAL 

81 

215 

.8 

6.0 

.6 

986.3 

26.4 

142 

STPAL 

81 

216 

.6 

7.4 

.  9 

984.3 

17.4 

142 

STPAL 

81 

217 

.8 

8.0 

2.8 

981.6 

15.0 

142 

STPAL 

81 

218 

.  5 

9.4 

4.0 

982.0 

16.0 

142 

STPAL 

81 

219 

.6 

6.3 

.7 

981.4 

11.8 

142 

STPAL 

81 

220 

.7 

5.3 

.2 

984.5 

10.6 

142 

STPAL 

81 

221 

1.0 

3.6 

.  5 

982.7 

16 . 0 

142 

STPAL 

81 

222 

1.0 

2.4 

-2.2 

975.5 

28.9 

142 

STPAL 

81 

223 

.7 

1.7 

-3.8 

973.1 

26.2 

142 

STPAL 

81 

224 

.  1 

2.4 

-5.8 

981 . 5 

14.1 

142 

STPAL 

81 

225 

.  1 

-.1 

-6  .  7 

991 . 5 

16.0 

142 

STPAL 

81 

226 

.6 

-1.6 

-6.5 

995.0 

20.1 

142 

STPAL 

81 

227 

1.0 

.6 

-2.4 

982.6 

18.1 

142 

STPAL 

81 

228 

1.0 

-.2 

-4.7 

984.8 

25.5 

142 

STPAL 

81 

3  1 

.7 

.4 

-6.3 

985.2 

23.2 

142 

STPAL 

81 

3  2 

.  1 

-5.1 

-12.9 
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18.3 

142 

STPAL 

81 

3  3 

.6 

-1.5 

-7.4 

980.7 

18.3 

142 

STPAL 

81 
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.3 

-1.3 

-7  . 1 

97  9.8 

16.7 

142 

STPAL 

81 

3  5 

.  3 

-1.9 

-10.1 

982.9 

16.2 

142 

STPAL 

81 

3  6 

.3 

-4.9 

-13.2 

98  9.4 

18.3 

142 

STPAL 

81 

3  7 

0.0 

-2.8 

-11.4 

997  .2 

8.1 

142 

STPAL 

81 

3  8 

.4 

.0 

-7.9 

998.3 

7.4 

142 

STPAL 

81 

3  9 

.8 

1.5 

-6 . 1 

996 . 7 

13.2 

142 

STPAL 

81 

310 

.7 

1.1 

-5.6 

995.0 

12.0 

142 

STPAL 

81 

311 

.0 

3.0 

-5.6 

993.8 

16.2 

142 

STPAL 

81 

312 

.2 

6.4 

-5.9 

984.6 

22.0 

142 

STPAL 

81 

313 

.3 

1.7 

-8.0 

988.1 

20.4 

142 

STPAL 

81 

314 

.  0 

4.7 

-7.4 

984 . 9 

22.9 

142 

STPAL 

81 

315 

.1 

7.6 

-4.7 

978.2 

27.1 

142 

STPAL 

81 

316 

.0 

2.9 

-7.3 

985.0 

16.0 

142 

STPAL 

81 

317 

.1 

1.9 

-10.3 

980.1 

25.0 

142 

STPAL 

81 

318 

.3 

-1.3 

-11.5 

981.1 

26.4 

142 

STPAL 

81 
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.6 

.8 

-9.2 

982.9 

25.9 

142 

STPAL 

81 

320 

.  0 

.8 

-9.2 

982.0 

13.2 

142 

STPAL 

81 

321 

.6 

3.2 

-7.2 

985.9 

10.4 

142 

STPAL 

81 

322 

.2 

5.9 

-5.6 

994.3 

8.1 

142 

STPAL 

81 

323 

.6 

6.7 

-5.3 

994.9 

7.4 

142 

STPAL 

81 

324 

.2 

8.3 

-5.4 

994.6 

13.4 

142 

STPAL 

81 

325 

.6 

6.7 

.1 

988.9 

14.1 

142 

STPAL 

81 

326 

.7 

5.4 

-.4 

987.4 

15.7 

142 

STPAL 

81 

327 

.8 

8.3 

-.6 

994.1 

23.8 

142 

STPAL 

81 

328 

1.0 

15.7 

6.5 

984.0 

25.9 

142 

STPAL 

81 

329 

1.0 

12.6 

9.2 

970.9 

17.4 

142 

STPAL 

81 

330 

1.0 

6.7 

2.6 

972.2 

13.9 

142 

STPAL 

81 

331 

.9 

8.7 

2.3 

971.9 

32.2 

:w2 

STPAL 

81 

4  1 

.7 

6.7 

.4 

978.7 

27.3 

142 

STPAL 

81 

4  2 

.5 

13.1 

2.0 

980.8 

17.4 

142 

STPAL 

81 

4  3 

.  9 

10.2 

4.5 

968.3 

24.5 

142 

STPAL 

81 

4  4 

1.0 

1.5 

-3.1 

971.5 

29.4 

U2 

STPAL 

81 

4  5 

.3 

2.4 

-8.3 

991.4 

12.7 

W2 

STPAL 

81 

4  6 

.  1 

7  . 1 

-2.7 

990.3 

26.2 

W2 

STPAL 

81 

4  7 

.  9 

11.9 

2.9 

983.8 

11.8 

W2 

STPAL 

81 

4  8 

.  9 

11.9 

-2.8 

985 . 7 

14.1 

W2 

STPAL 

81 

4  9 

.5 

11.6 

-1.3 

989.9 

20.8 

W2 

STPAL 

81 

410 

.5 

15.5 

2.5 

987.5 

17.6 

W2 

STPAL 

81 

411 

.  9 

8.7 

-1.9 

991.4 

18.5 

W2 

STPAL 

81 

412 

.7 

11.6 

4 . 9 

990.8 

19.7 

W2 

STPAL 

81 

413 

1.0 

10.4 

6.5 

990.2 

20.6 

W2 

STPAL 

81 

414 

0  .  0 

3.8 

-7.4 

10  05.4 

24.1 

W2 

STPAL 

81 

415 

.  0 

8.4 

-6.3 

1006.5 

21.8 

U2 

STPAL 

81 

416 

.8 

10.9 

3.5 

990.3 

20.6 

142 

STPAL 

81 

417 

.4 

13.6 

1.4 

985.9 

23.8 

142 

STPAL 

81 

418 

.  5 

8.3 

-1.0 

996.7 

12.7 

142 

STPAL 

81 

419 

.8 

10.4 

3.8 

995.3 

18.8 

142 

STPAL 

81 

420 

.4 

5.3 

-6 . 3 

999.9 

17.6 

142 

STPAL 

81 

421 

1.0 

5.6 

-1 . 3 

987 . 1 

17.8 

142 

STPAL 

81 

422 

1  .  0 

11.0 

8.1 

978 . 0 

11.8 

142 

STPAL 

81 

423 

1  .  0 

5.1 

1.8 

976.6 

30.3 

142 

STPAL 

81 

424 

.7 

7.8 

.  6 

985.3 

18.8 

142 

STPAL 

81 

425 

.8 

9.2 

1.9 

985.3 

10.6 

U2 

STPAL 

81 

426 

.6 

14.1 

4.9 

987 . 4 

10.2 

142 

STPAL 

81 

427 

.8 

16.5 

9 . 6 

986.5 

21.8 

142 

STPAL 

81 

428 

.8 

10.1 

4.9 

985.5 

17  .4 

142 

STPAL 

81 

429 

.  9 

9.2 

5.3 

982.8 

10.2 

142 

STPAL 

81 

430 

.4 

9  .  3 

5.2 

987 . 0 

10.4 

142 

STPAL 

81 

5  1 

.4 

9.0 

-1.5 

992.8 

12.7 

142 

STPAL 

81 

5  2 

.4 

11.6 

2 . 3 

986.3 

16.4 

142 

STPAL 

81 

5  3 

1.0 

14.3 

11.0 

978.9 

18.5 

142 

STPAL 

81 

5  4 

1 .  0 

15.1 

12.4 

983.8 

15 . 5 

142 

STPAL 

81 

5  5 

.6 

12.8 

4.3 

995.8 

20 . 1 

142 

STPAL 

81 

5  6 

.3 

11.1 

-1.0 

1000.2 

8.8 

142 

STPAL 

81 

5  7 

.7 

11.5 

.9 

996.2 

6.3 

142 

STPAL 

81 

5  8 

.8 

12.3 

3.3 

987.7 

13.4 

142 

STPAL 

81 

5  9 

.8 

8.8 

2.3 

991 . 1 

20.6 

142 

STPAL 

81 

510 

.0 

6.9 

-5.7 

995.0 

21.3 

142 

STPAL 

81 

511 

.  1 

9.5 

-6.1 

993.2 

16.9 

142 

STPAL 

81 

512 

.8 

10.5 

-4 . 1 

992.6 

14.1 

142 

STPAL 

81 

513 

.2 

13.3 

-.0 

992.6 

10.2 

142 

STPAL 

81 

514 

.1 

14.4 

.2 

989.7 

6.3 

142 

STPAL 

81 

515 

.3 

16.3 

1.3 

989.2 

6.0 

142 

STPAL 

81 

516 

.9 

17.3 

3.5 

992.9 

14.6 

142 

STPAL 

81 

517 

1.0 

13.3 

4.8 

997.7 

19.4 

142 

STPAL 

81 

518 

.1 

13.1 

-3.0 

999.0 

18.5 

142 

STPAL 

81 

519 

0.0 

14.4 

-1.3 

997.1 

7.6 

142 

STPAL 

81 

520 

0.0 

17.7 

-.1 

995.2 

8.8 

142 

STPAL 

81 

521 

.4 

20.0 

1.8 

991.3 

19.9 

142 

STPAL 

81 

522 

.8 

21 . 0 

10.4 

985.3 

26.9 

'142 

STPAL 

81 

523 

1.0 

18.9 

14.2 

981.6 

17.4 

142 

STPAL 

81 

524 

1.0 

13.5 

9.7 

980.6 

18.3 

142 

STPAL 

81 

525 

1.0 

12.6 

9.2 

983.2 

14.4 

142 

STPAL 

81 

526 

1.0 

15.3 

11.6 

988.8 

10.2 

142 

STPAL 

81 

527 

1.0 

16.6 

10.5 

990.6 

8.8 
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Table  B1 

Initial  Conditions  Used  for  1981  and  1982  Eau  Galle  Simulations 


Variable 


Units 


Model 

Acronym 


Algae  1 

g/m3 

ALGAE (  ,1) 

0 

0 

Algae  2 

g/m3 

ALGAE (  ,2) 

0 

0 

Algae  3 

g/m3 

ALGAE (  ,3) 

2.77 

0.33-0.38 

Temperature 

°C 

TEMP 

5. 8-6.0 

6.9-7.  1 

Alkalinity 

g/tn3 

ALKA 

108.0-110.0 

80.0-82.0 

Ammonia-N 

g/m3 

CNH3 

0.033-0.052 

0.  184  -0.200 

Nitrite  nitrate-N 

g/m3 

CN02 

0.654-0.682 

1.14-1.17 

Refractory  dissolved 
organics 

g/m3 

RFR 

13.5-15.3 

10.6-12.6 

Labile  dissolved 
organics 

g/m3 

DOM 

5.79-6.54 

4. 5-5.4 

Detritus 

g/m3 

DETUS 

2.4 

2.4 

Oxygen 

g/m3 

OXY 

12.7-13.3 

10.3-10.4 

Orthophosphate-P 

g/m3 

P04 

0.110 

0.081-0.083 

Total  dissolved  solids 

g/m3 

TDS 

157.0-167.0 

133.0-153.0 

Zooplankton 

g/m3 

ZOO 

0.1 

0.1 

pH 

DL* 

PH 

8. 0-8. 2 

7.8 

Suspended  solids 

g/m3 

SSOL 

23.4-24.6 

11.3-13.3 

Sediment 

g/m2 

SEDMT 

501.0-5001.0 

501.0-5001.0 

Fish 

kg/ha 

FISH 

55.0 

55.0 

Oxidized  manganese 

g/m3 

CMN4 

0.1 

0 

Reduced  manganese 

g/m3 

CMN2 

0.1 

0 

Oxidized  iron 

g/m3 

FE3 

0.8-1. 2 

0.5 

Reduced  iron 

g/m3 

FE2 

0.1 

0 

Iron  sulfide 

g/m3 

FESB 

0 

0 

Sulfate 

g/m3 

S04 

7. 1-9.3 

11.1-12.9 

Sulfide 

g/m3 

S2 

0 

0 

dimensionless . 
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Table  B 1  (Cone luded) 


Model 

Variable _  Units  Acronym  1  98 1 _  _ 1  982 


Sediment  iron 

g/m3 

FE 

1310.0 

1310 

Sediment  iron  sulfide 

g/m3 

FESA 

0 

0 

Sediment  sulfur 

g/m3 

S 

1.6 

1.6 

Sediment 

orthophosphate-P 

g/tn3 

XP04 

102.0 

102.0 

Sediment  nitrogen 

g/m3 

CN 

1210.0 

1210.0 

Silica 

g/m3 

SI 

5.52-6.10 

0.33-0.38 

Colif orms 

e/m3 

COLIF 

220.0-300.0 

111.0 

1982  Eau  Galle  Simulations 


Table  Cl  (Continued) 


production  rate 


Table  Cl  (Continued) 
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Table  CL  (Continued) 


Model  Acron 


Table  Cl  (Continued) 
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Anerobic  (Continued 


